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Abstract 

Using continued fractions we study the ground state properties of the spin-| Ising chain in a transverse 
field with periodically varying interaction strengths and external fields. We consider in detail the chain having 
the period of modulation of interactions equals 2 and compare the results obtained with those corresponding 
to the spin-i isotropic XY chain in a transverse field. In contrast to the behaviour of the transverse XY 
chain, the transverse Ising chain does not exhibit a step-like magnetization vs. field dependence caused 
by the alternation of bonds, its susceptibility exhibits a logarithmic singularity at the field determined by 
interaction strengths, and it is stable with respect to spin-Peierls dimerization. 
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Spin 



i XY chains provide an excellent ground for a rigorous study of different properties of the low- 



dimensional quantum magnetic systems due to the fact that with the help of the Jordan- Wigner transformation 
such spin models can be mapped onto noninteracting spinless fermions and as a result many statistical mechanics 
calculations can be performed exactly. 1 In what follows we shall consider the periodic alternating XY chain in 
a transverse field with extremely anisotropic spin coupling (in short called here the transverse Ising chain) in 
order to study the generic effects induced by regular alternation of the nearest neighbour interaction strengths 
and the external fields. The periodic alternating chain can be viewed as a chain having several sublattices. 
There are few papers dealing with the statistical mechanics properties of XY chains on two sublattices. 2-5 
We want to emphasize that in the above mentioned limit of the transverse Ising chain the calculation of the 
thermodynamic quantities can be performed in a rather general manner covering all chains with an arbitrary 
period of alternation. The present paper can be viewed as a nontrivial extension of the recent study of the 
thermodynamics of the spin-^ isotropic XY chain in a transverse field (called here the transverse XY chain) 
with regularly alternating fields and bonds. 6 In particular, we compare below the zero temperature dependences 
transverse magnetization vs. transverse field and static transverse susceptibility vs. transverse field for the 
regularly alternating transverse Ising and XY chains. On the other hand, it is known that the transverse XY 
chain is a simple system exhibiting a spin-Peierls instability. An analysis of the ground state energy of the 
periodic alternating transverse Ising chain of period 2 allows one to examine a spin-Peierls instability of that 
chain with respect to dimerization. We find that the transverse Ising chain is stable with respect to dimerization 
which demonstrates a role of the anisotropy in spin coupling for a spin-Peierls instability. 

We consider the N — > oo spins ^ on a ring with the Hamiltonian of the nonuniform transverse Ising model 



We assume regular nonuniformity in (Q), i.e. the transverse field Q n at the site n as well as the exchange coupling 
I n between the neighbouring sites n and n + 1 vary regularly from site to site with period p, i.e. we have a 
sequence of parameters fiiii 02-^2 ■ ■ ■ QpIpQiIiQzh ■ ■ ■ Qplp ■ ■ ■ ■ 

Our goal is to examine the thermodynamic quantities of spin model (Q). In order to do this, we shall first 
express the Hamiltonian in terms of fermion operators. This can be done in the usual way by applying the 
Jordan- Wigner transformation. 1 As a result one gets a model of spinless fermions on a ring described by the 
Hamiltonian 

N ( 1\ 1 N 

H = ^ n « f C+C„ - - J + - Y In (c+ C n+l + C n C n+l ~ C n C+ +l - C„C, i+ i) , C+ +1 = cf , C N+1 = C X (2) 
n=l > ' n=l 

(the boundary term, that is unimportant as far as the thermodynamics is concerned, has been omitted). Unlike, 
in Ref. 6, we cannot proceed directly by using a continued fraction representation for the diagonal one-fermion 
Green functions of the tight-binding spinless fermions since the Hamiltonian (0) contains the products of two 
creation (annihilation) operators. However, it is well known (sec, e.g., 1 ' 11 > 12 ) that after introducing new 




(1) 



operators % = J2iLi (9kiCi + h kl cf), 77^ = J2iLi ( h kiCi + gkicf) the Hamiltonian (§) transforms into 
N ( 1\ 

H = Xj A fc [vim - 2 J . {vt ,Vg} = Skg, {Vk,Vq} = l^/fe = (3) 
if the unknown coefficients = \{<&ki + ^fei); ^fej = hi^ki ~ ^h) are determined from the following equations 

o n _ii n _ 1 $ fe) „_ 1 + (n£ + _ A 2) $fen + ^ n / n $ fc n+1 = o ; $ fc0 = $ few , $ fciAr+1 = $ fel ; 

nn/n-l*M-l + + J n - A fe) *fen + ^n+l^n^fe.n+l = 0, f kO = ^ kN , ^k,N+l = *fcl< (4) 

Eqs. (Q) formally coincide with those describing displacements of particles in a nonuniform harmonic chain 
with nearest neighbour interactions. To find the distribution of the squares of "phonon" (magnetic excitation) 
frequencies R{E 2 ) = i J2k=i 5 { E% - A D 

we may use the Green function approach. Consider, for example, the 
Green functions G nm = G nm (E 2 ) that satisfy the set of equations 

( E — l) ^ nm ^n— X^n— l^"n— l,m ^n^ti^nll.m — $nm (5) 

with periodic boundary conditions implied. Knowing G nm one immediately finds the density of states R(E 2 ) 
through the relation 

JV 

,T.V 



1 N 

R{E 2 ) =T-r F ^ImG„„(^ 2 ±ze), e -> +0. (6) 



n=l 

r"2\ 



As follows from (|3|) all the thermodynamic quantities can be expressed through R(E ). For example, the 
Helmholtz free energy per site is given by 



2 r°° / 
f = -- dEER(E 2 )\n 2 cosh 
P Jo \ 



f3E 



(7) 



Obviously, R(E 2 ) can be obtained with the help of the Green functions introduced on the basis of the set of 
equations for the coefficients ~$>kn (||)- We have performed such a calculation and found the same result for 
R{E 2 ) in the cases considered below (Eqs. (||) and (|ll|)). 

We have now to calculate the diagonal Green functions G nn ■ Let us use the continued fraction representation 
for G nn that follows from (||) 

1 



Gn, n — 



Z?2 _ 02 _ j2 _ A _ A + 

O 2 T 2 



e 2 - n 2 



n-1 A n-2 E 2 -il.' 2 „-/ 2 ,- 



O 2 J 2 

E 2 - Q, 2 - II " +1 " +1 



n 2 i 2 

nn -. (8) 



n+l ± n E 2 -nl +2 -I 2 + 1 - 



For any finite period of varying f2 n and /„ the continued fractions in (Bh become periodic (evidently, the limit 
N — > oo is hinted) and can be easily calculated by solving quadratic equations. As a result we get rigorous results 
for the Green functions, density of states and thermodynamic quantities (|7|) of the periodic alternating spin 



model ([y). Note that the thermodynamic quantities are not sensitive to the signs of f2 n and 7„. Therefore, one 
may assume £l n > 0, 7„ > without loss of generality. 

It should be remarked at this point that the possibility of obtaining R(E 2 ) exactly in the manner described 
above exists only for the transverse Ising chain. For the transverse XY chain with an arbitrary anisotropic 
exchange coupling one arrives at a set of five diagonal coupled equations (not three, as above), corresponding 
to a nonuniform harmonic chain with nearest and next nearest neighbour interactions. 1 ' 12 Hence, one cannot 
proceed as above. 

To illustrate how the method works we begin with the uniform case, namely, the periodic alternating chain 
having period 1. For such a case, 

1 



A" = A+ = -(E 2 - n 2 - I 2 ± sJ(E 2 - n 2 - I 2 ) 2 - 4Q 2 I 2 ), 
G 1 



y/(E 2 -n 2 -i 2 ) 2 -4n 2 i 2 



and therefore 



=, if ai < E 2 < a 2 , 



R(E- J )={ * y/-^- ai )(E'-a 2 y (9) 

0, otherwise, 
with ai = (n — I) 2 , a,2 = (n + I) 2 - From Eq. (||), the ground state energy per site, e = — / °° dEE 2 R(E 2 ), is 

1 dEE 2 

eo = - 



W^r y/-(E 2 ~a 1 )(E 2 -a 2 )' 



which after the change of variable E = y a 2 — (a 2 — a>i ) sin 2 0, leads to the known 13 expression 

e o 1 P u 17, , xn2 77 . 2 ] x ^ 



7 7T 



(1 + A) -4Asin 2 0, A=-. (10) 



Let us now turn to the periodic alternating chain having period 2. Following the procedure above, one finds 



1 [2g^— 0£-0|— j£j . f n(T?2\ ^ n 



(11) 



0, otherwise, 

b(e 2 ) = 4n 1 fi 2 ,7 1 2 7 2 2 - (e a - (n'l + n 2 2 + if + 1\) e 2 + njnl + ifif) 2 

= -(E 2 - b^E 2 - b 2 )(E 2 - b 3 )(E 2 - 64), 

{M - [\ (nf + ni + A 2 + ^ 2 ± 7 (n? + nl + 7 2 + 7 2 ) 2 - 4(0^2 ± hh) 2 ^ I . 

Obviously, ([ll]) recovers the result for the uniform chain (^J) if Qi = fia = ^, 7i = 72 = 7. The resulting 
density of states R(E 2 ) ( |ll|) yields all thermodynamic quantities of the regularly alternating transverse Ising 
chain r2 1 7 1 r2 2 72^2i7 1 fl 2 72 • • • • 

Consider further, in some detail, the alternating transverse Ising chain with p = 2 with the dimerization 
ansatz Qi = O2 = 7i = 7(1 + S), I2 = 7(1 — 6) in mind, where < 6 < 1 is the dimerization parameter. For 
such a chain the ground state energy according to (|Tl| ) can be written as 

- 1 ■ Jj dej) Ufl + A 2 + 5 2 + 2VW+^J l + \ 2 +S 2 - 2Vw^j , 

W = X 2 + 5 2 - \ 2 (l-6 2 )sm 2 (f>. (12) 
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If 5 = (|l2| ) transforms into the ground state energy of the uniform chain (|T 

Further, we calculate the zero temperature transverse magnetization m z = ^ by differentiating the r.h.s of 
Eq. ( p"2| ) with respect to A. In Fig. 1 the obtained dependence of the transverse magnetization on the transverse 
field for several values of the dimerization parameter S is plotted. As one can see from Fig. 1 the transverse 
Ising chain with regularly alternating bonds does not exhibit a step-like dependence of magnetization vs. field 
and the magnetization profile only smoothly deforms with the increasing 5. This behaviour is contrary to that 
of the transverse XY chain (for which the Hamiltonian (Q) would contain s^s^ +1 + s^s^+i instead of s^sjj +1 ), 
since the latter model does show a step-like dependence of magnetization vs. field (i.e. a plateau at m z = 0) 
for regularly alternating bonds 6 (the light curves in Fig. 1). The dissimilarity is conditioned by the different 
symmetries of these models, since Y) n — i s z n does not commute with the Hamiltonian of the transverse Ising 
chain but it does commute with the Hamiltonian of the transverse XY chain. It is worth to consider in addition 
both models in the limit 6 = 1 when the chain splits into noninteracting clusters consisting of two sites. The 
zero temperature transverse magnetization follows from the formula m z = i(G>S|sf + sf IGS 1 ) where \GS) is the 
ground state eigenvector of the two-site cluster Hamiltonian. For the transverse Ising model \GS) smoothly 
varies with increasing fi. Moreover, (GS\sf + s%\GS) = if Q = and (GS\sf + s^\GS) -> -1 if fi -> oo. 
Unlike, for the transverse XY model the ground state is singlet (and thus (GS\sf + s^GS) = 0) if < ft < 21. 
However, if f2 exceeds 21 the ground state becomes triplet and (GS\sf + s^GS 1 ) abruptly changes to —1. 

The ground state static transverse susceptibility \ zz = ^f- = j J^? e "j l5 - > follows straightforwardly from Eq. 
(|l2|). In the left panel in Fig. 2 the obtained dependence of the static transverse susceptibility on the transverse 
field for several 6 is plotted. Using Eq. (|lj) it is easy to show that \zz for the dimerized chain consists of a 
finite part and a singular part 

A 2 H (1- (l-5 2 )sin 2 ,' N2 



\ zz = finite term / dip — - — . (13) 

27r Jo wi\/l + y + 6 2 -2VW 

As it follows from (|l^) the ground state static transverse susceptibility exhibits a logarithmic singularity at 



A = A* = Vl — S 2 , i.e. Xzz ~ hi |A — A*|, A — > A*. It should be noted, that \zz for the corresponding dimerized 
XY chain exhibits two square root singularities at A^ = 28 and A£ = 2 (right panel in Fig. 2, see also Ref. 6). 

To discuss a spin-Peierls dimerization in the adiabatic limit one must examine the dependence of the total 
energy £ [S) = eo(S) + aS 2 , which consists of the magnetic part eo(S) jl2|) and the elastic part aS 2 , a > 0, on the 
dimerization parameter 5. In Fig. 3 one can see the dependence of the magnetic energy ( |l2| ) on 5 for different 
values of the transverse field. From these plots one concludes that the magnetic energy as a rule decreases with 
the increase of 5 (except for A = and A — > oo when ^j^- (|l^) does not depend on 5 and equals —\ and — -|, 
respectively) that is a necessary condition for the existence of the dimerized phase. However, the magnetic 
energy decreases too slowly in comparison with the increase of the elastic energy that provides stability of the 
chain with respect to dimerization. To demonstrate this, we calculate a = —-ggz eoi f^ by differentiating the 
r.h.s. of Eq. ([l2|) with respect to S 2 , and display the dependence of this quantity on S for several values of the 
transverse field in Fig. 4. These plots show that £(S) may exhibit a maximum (but not a minimum) at the 



nonzero value of dimerization parameter for lattices having small j as the slope of the curves a = — e ") 
vs. 6 manifests. Hence, the uniform chain with <5 = is favourable. On the contrary, for the XY chain the 
ground state energy decreases sufficiently rapidly to provide a stability of the dimerized phase. In Fig. 4 we 
also display the dependence of a — —-J^i^p- on <5 for the transverse XY chain (light curves) (see also Ref. 6) 
emphasizing explicitly the difference between those two spin models. 

It is worthy to note, that the absence of a spin-Peierls dimerization in the transverse Ising chain may be 
anticipated on the basis of the following reasoning. 14 Considering the transverse XY chain one notes that 
for f2 = the ground state corresponds to a half-filled fcrmion band, the system is gapless and the Peierls 
mechanism does work. As f2 increasing becomes equal to / the ground state corresponds to an empty fermion 
band, the system becomes gaped and the dimerized phase does not appear. The ground state of the transverse 
Ising chain 13 always corresponds to an empty fermion band, the system is gapped unless £1 = / and as a result 
no Peierls dimerization should be expected. 

To conclude, we have proposed to use continued fractions for the calculation of thermodynamic quantities 
of regularly alternating spin-i transverse Ising chains. For such a chain, in which the bonds are alternating 
in magnitude along the chain with period 2, we have found no plateaus in the dependence magnetization vs. 
field. The susceptibility of the chain exhibits a logarithmic singularity at the field determined by interaction 
strengths. Moreover, we have found that in the adiabatic limit the transverse Ising chain is stable with respect to 
a spin-Peierls dimerization. These findings differ dramatically from the corresponding results for the transverse 
isotropic XY chain, and therefore it would be of great interest to study the periodic alternating transverse 
anisotropic XY chain that connects both these limiting cases. 
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List of figure captions 

FIGURE 1. The ground state transverse magnetization vs. transverse field for the transverse Ising chain with 
regularly alternating bonds with the period 2. 5 = (uniform chain), 6 = 0.3, 6 = 0.6, and 5 — 0.9 correspond 
to solid, long-dashed, short-dashed, and dotted bold curves, respectively. For comparison the corresponding 
results for the transverse XY chain are shown by light curves. 

FIGURE 2. The ground state static transverse susceptibility vs. transverse field for the transverse Ising 
chain (left panel) and the transverse XY chain (right panel) with regularly alternating bonds with the period 
2. S — (uniform chain), 5 = 0.3, 5 = 0.6, and 5 = 0.9 correspond to solid, long-dashed, short-dashed, and 
dotted curves, respectively. 

FIGURE 3. The change of the ground state magnetic energy d = — Sei2i vs . ,5 f or the transverse Ising 
chain for A = 0, 0.25, 0.5, 1, and 2 (solid, long-dashed, short-dashed, dotted, and dashed-dotted bold curves, 
respectively). For comparison the corresponding results for the transverse XY chain are shown by light curves. 

FIGURE 4. The dependence a = —^i eo ^p vs. 8 for the transverse Ising chain with A = 0, 0.25, 0.5, 1, 
and 2 (solid, long-dashed, short-dashed, dotted, and dashed-dotted bold curves, respectively). For comparison 
the corresponding results for the transverse XY chain are shown by light curves. 
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Figure 1: The ground state transverse magnetization vs. transverse field for the transverse Ising chain with regularly 
alternating bonds with the period 2. S = (uniform chain), S = 0.3, S = 0.6, and S = 0.9 correspond to solid, long- 
dashed, short-dashed, and dotted bold curves, respectively. For comparison the corresponding results for the transverse 
XY chain are shown by light curves. 




Figure 2: The ground state static transverse susceptibility vs. transverse field for the transverse Ising chain (left panel) 
and the transverse XY chain (right panel) with regularly alternating bonds with the period 2. 8 = (uniform chain), 
S — 0.3, 8 — 0.6, and 8 — 0.9 correspond to solid, long-dashed, short-dashed, and dotted curves, respectively. 



d 




Figure 3: The change of the ground state magnetic energy d — — SaLi vs g f or the transverse Ising chain for 

A = 0, 0.25, 0.5, 1, and 2 (solid, long-dashed, short-dashed, dotted, and dashed-dotted bold curves, respectively). For 
comparison the corresponding results for the transverse XY chain are shown by light curves. 
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Figure 4: The dependence a 
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vs. S for the transverse Ising chain with A = 0, 0.25, 0.5, 1, and 2 (solid, 
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long-dashed, short-dashed, dotted, and dashed-dotted bold curves, respectively). For comparison the corresponding 
results for the transverse XY chain are shown by light curves. 



